Quasienergy spectra of a charged particle in planar honeycomb lattices 
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The low energy spectrum of a particle in planar honeycomb lattices is conical, which leads to 
the unusual electronic properties of graphene. In this letter we calculate the quasienergy spectra 
of a charged particle in honeycomb lattices driven by a strong AC field, which is of fundamental 
importance for its time-dependent dynamics. We find that depending on the amplitude, direction 
and frequency of external field, many interesting phenomena may occur, including band collapse, 
renormalization of velocity of "light", gap opening etc.. Under suitable conditions, with increasing 
the magnitude of the AC field, a series of phase transitions from gapless phases to gapped phases 
appear alternatively. At the same time, the Dirac points may disappear or change to a line. We 
suggest possible realization of the system in Honeycomb optical lattices. 
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Lattice structure has important impact on the dynam- 
ics of a system through the related band structure. For 
many materials with square/cubic lattice, their energy 
spectra near the bottom of conduction band are parabolic 
and the low energy excitation is a quasiparticle with an 
effective mass. However graphene, a two-dimensional 
monolayer of carbon atoms, has a linear dispersion rela- 
tion near so-called Dirac points due to its special honey- 
comb lattice[l|. Therefore in the long distance limit, the 
basic excitation is a massless Dirac Fermion. There are 
extensive ongoing experimental and theoretical studies 
on graphene On one hand graphene provides an ex- 
cellent condensed-matter analogue of (2-1-1) dimensional 
quantum electrodynamics, where the effective velocity of 
"light" is 1/300 of the real velocity of light. On the 
other hand the unusual energy spectrum also leads to 
novel transport properties, such as anomalous quantum 
Hall effect [3j]. Besides graphene, the optical honeycomb 
lattice has also been realized in experiments The 
flexible tunability of optical lattices gives people much 
more opportunities to explore the interesting physics in 
two-dimensional systems with honeycomb lattices. Re- 
cently much attention has been drawn to this direction, 
for instance Haldane's quantum Hall effects was proposed 
in ultracold atoms in optical honeycomb lattice [5|, P- 
orbital physics in the honeycomb optical lattice was also 
addressed Bose-Einstein condensation in a honey- 
comb optical lattice was considered in [7], just name a 
few. 

An important issue here is the adjustment of band 
structure. One convenient and effective way is through 
external time-dependent field. For a system driven by 
a time-periodic electric field, band suppression Q and 
band collapse[l3| appear. There are also many other 



interesting phenomena, including coherent destruction 
of tunneling fs'], dynamical localizationfll'l, photovoltalic 
effects jl 21 ■ and photon- assistant Fano resonance 13| etc.. 
Many of them have been observed in electronic system 
and/or optical lattices. As the energy spectrum is es- 
sential for studying the optical and electronic properties 
of the relevant system, the quasienergy spectrum plays a 
key role in understanding the time-dependent phenom- 
ena. 

Thus it is of great importance and interest to study the 
quasienergy spectra of a particle in honeycomb lattice in 
the presence of an AC field. This issue is addressed in 
detail in this paper. We show that the quasienergy spec- 
trum of honeycomb lattice has a quite rich structure. A 
series of phase transitions (from gapless phase to gapped 
phase) appear with changing of the amplitude, direction, 
and frequency of the external electric field, which is quite 
easy in experiments. In particular, under special condi- 
tions the particle can be localized in one direction, or 
two directions, thus the dimension of the system reduces 
effectively from two dimension (2D) to one dimension 
(ID) and to zero dimension (OD). The changing of exter- 
nal field may also lead to the effective mass generation 
and renormalization of the velocity of "light" . 

The planar honeycomb lattice (Fig. 1(a)) consists of 
two equivalent sublattices A (black dots) and B (white 
dots). We consider the charged particle hopping between 
the nearest neighbor sites from different sublattices. In 
addition, the system under consideration is subject to a 
time-periodic electric field E(<) = Eocos(w<). The tight- 
binding Hamiltonian has the form 



eE(t) 



(1) 



*Email address: zhang _wei@iapcm.ac.cn 



where < i,j > refers to the nearest neighbor lattice sites, 
Tj is the position of site i. The probability amplitude 
a{ri) for a particle at site i of sublattice A (or /3{ri) at 
site of sublattice B ) can be calculated by solving the 



FIG. 1: (a) Planar honeycomb lattice with electric field along 
Y direction. Black dots are sites of sublattice A and white 
dots are sites of sublattice B. (b) The effective dimerization 
bond pattern. The effective hopping constants along the (red) 
thick bonds are different from those along (black) thin bonds. 
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FIG. 2: Red solid curve; |^1/2|B1 = \Jo(x) /2Jq{x /2)\ versus 
x; Blue dotted curve; |i/2Jo(a;)| versus x. 



corresponding Schrodinger equation. It is very helpful to 
make the unitary transformation a(rj) a(ri)e~*^*^*-' ''' 

(/3(r,) ^ ^(r,)e*^(*)-'--), where A{t) = ^^sin{ujt), 
with a the lattice constant. The subsequent Fourier 
transformation leads to 



the quasienergy spectra around the minimum of conduc- 
tion band/maximun of valence band for different values 
of eEoa/uj, which are obtained by exact numerical calcu- 
lation based on Eqs. (2) and (3). 
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Due to the Floquet theorem, the quasienergy can be 
obtained by diagonalizing the time evolution operator 
T dte~'^^^^'*\ where T refers to time ordering and T 
is the period of the AC field . We first consider the case 
of electric field along Y direction. In the high frequency 
limit V/io << 1, we can obtain the following analytical 
result for quasienergy spectrum 




s = ±V\ 



AB^{cos—k^a)^ 



AAB cos ——Kxa cos - 
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FIG. 3; Quasienergy spectra for eEoa/ij—3.5 (in regime I), 
eEoa/Lj—4.60 (in regime II), ei?offl/w=4.81 (band collapse 
^condition), eEoa/ui=5.52 (in regime III), ei5oa/t<;=10.8 (in 
k a'PSY"^^ eEoa/u)=12.5 (in regime V), are obtained by 

^ exact numerical calculation. V/uj = 0.02. 



where A = Jo{eEoa/uj) and B = Jo{eEoa/2u;), Jq is the 
zero-th order Bessel function. When Eq ^ 0, A, B ^ 1, 
we reproduce the results for energy spectrum of a system 
with honeycomb lattice without AC electric field. 

Quantum phase transition The conduction band 
and valence band may form conically shaped valleys that 
touch at some specific points (called conical points or 
Dirac points). These Dirac points are determined as 

kyU = 27r/3, cos(^^fc2;a) = A/2B. It is clear that there 
is no solution, when \A\ > 2\B\. In this situation, an 
energy gap is generated. In Fig. 2 (red solid curve) we 
show |^|/2|B| = I Jo(x)/2Jo(a:/2)| versus x = eEoa/u. 
From Fig. 2 it can be seen that gap opening happens in 
the regimes (for eEoa/u) II=[4.1, 5.1], IV=[10.3, 11.3],.... 
Figures 3(a)-(f) (corresponding to regimes I to V) show 



We can see from figure 3 that as eEoa/uj increases, 
the system undergoes a series of phases corresponding to 
the regimes I-V, i.e., gapless phase (regime I) gapped 
phase (regime II) — > gapless phase (regime III) — > gapped 
phase (regime IV) gapless phase (regime V). The gap 
opening introduces a new energy scale and may also 
lead to the generation of an effective mass as seen in 
Figs. 3(b) and (e), since the quasienergy spectrum near 
the bottom of conduction band is parabohc. We find 
the gap and the effective mass are proportional to D = 
{\A\-2\B\)e{\A\ - 2\B\). Thus D is the order parameter, 
which is zero in the massless phases and nonzero in the 
gapped phases. As one can see from figure 3, with tuning 
AC field, the Dirac points may disappear in the gapped 
phase, or change from discrete points to a line as shown 
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in Fig. 3(d). Near e = 0, the density of quasienergy 
states D{e) is zero for gapped phases, D{e) ~ \e\ near 
discrete Dirac points, and D{e) ^ constant near the hne. 
According to the generahzcd Kubo-Greenwood formula 
for time-dependent system[14)], the conductivity depends 
on the density of quasienergy states. Therefore by simply 
tuning the strength of external AC field, we can realize 
the transitions from "metal" phase to "semimetal" phase 
and to "semiconductor" / "insulator" phase. 

Band collapse A typical character of Fig. 2 is that 
|^|/2|B| is divergent at eE^ajuj = 4.81,11.0,.... By 
checking Eq. (4), it is easy to see that band collapse 
happens and a gap opens under this condition B=0 or 
Jo{eEoa/2u}) = 0. Figure 3(c) shows the quasienergy 
spectrum for eEoa/oj = 4.81, twice of 2.405, the first 
zero point of Jq. As expected, we see the flat conduction 
and valence bands, which are drastically different from 
the remaining figures of Fig. 3, and the case without an 
AC field. Moreover we see a gap between conduction and 
valence bands opens. From eq. (4), one can also see that 
the quasienergy band collapses in the Y direction, when 
A=0, i.e., eEoa/uj = 2.405,.... 

Time-dependent dynamics The collapse of 
quasienergy spectrum implies dynamical localization. To 
verify this picture, we calculate the time evolution of a 
particle initially at the origin. In our calculation, the sys- 
tem size is 101 x 101, time stcp=0.02 (in unit of 1/w), 
and the hopping constant is V/uj = 0.02 as that in the 
calculation of quasienergy. In the long time regime, the 
mean square displacement of the particle shows the time- 
dependent behavior as t". In figure 4(a) we show the 
time dependence of the mean square displacement. We 
can see that when the system is away from band collapse 
condition {eEoa/uj = 4.81), a = 2, the dynamics is a 
typical ballistic motion. Under the band collapse con- 
dition, the mean quare displacement is finite (< 1) for 
all the time, i.e., a = 0, indicating the dynamical lo- 
calization of the particle due to the presence of an AC 
field. In Fig. 4(b), we show the time dependence of the 
mean square displacement in the X and Y directions for 
eEfja/uj — 2.405, under which the quasienergy spectrum 
is flat in the Y direction. It is clear that the particle is 
localized in the Y direction (< >< 0.18), indicating 
a dimension reduction from 2D to ID. 

Renormalization of velocity of "light" For sys- 
tem in the gapless phases, we may obtain the velocity of 
"light" by expanding the quasienergy spectrum around 
Dirac points 



(5) 



where {kxO,, kyu} refer to the Dirac points. We can see 
that the velocity of "light" can be tuned by the amplitude 
and frequency of external electric field. From Fig. 5(a), 
one can see that the velocity of "light" for eE^a/oj = 3.5 
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FIG. 4: (a) Time dependence of the mean square- 
displacement for eEoa/u) = 4.81 and eEoa/u = 5.50; (b)Time 
dependence of the mean square-displacement in the X and Y 
directions for eEoa/u) — 2.405. 
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of that without external electric field. 



The external field may also lead to the anisotropy of ve- 
locity of "light", see figure 5(b) and Eq. (5). An ex- 
tremal example is that the velocity in one direction is 
finite, while in the other direction becomes zero under 

the condition A=0, i.e. cos(^^fcj;a) — 0, as shown in 
Fig. 3(d). The tunability of velocity of "light" gives us 
more opportunities to explore relativistic effects in elec- 
tronic systems or optical lattices. 
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FIG. 5: Renormalization of the velocity of "light", 
(a) Quasienergy spectrum for eEoa/u>—0.0 and 3.5. (b) 
Anisotropy of velocity of "light". Quasienergy spectrum for 
eEoa/uj = 12.5. For the purpose of comparison, we have 
shifted the coordinates of Dirac points. 

Dependence on the direction of electric fields 

The quasienergy spectrum is dependent on the direc- 
tion of electric field. For the case of an electric field 
along X direction, the quasienergy spectrum still takes 

the form of Eq. (4) with A=l and B = Jo{^eEoa/u!). 
In this case the condition for the existence of Dirac points 

is \2jQ{^eEQa/uj)\ ^ 1. From Figure 2(b) (dotted 

blue curve), we see that when X^eEoa/uj > 1.52 or 
eEoa/uj > 1.76, a gap always opens, which is quite con- 
trast to the case with electric field along Y direction dis- 
cussed above. The band collapse appears under the new 

AT 

condition Jo{X£-eEoa/uj) — 0. 

In the more general situation, the quasienergy spec- 
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trum has the form £ — ±t^^| aje^*"'' "^ P, where 

aj = Jo(eaEo ■ Uj/w). The general condition for the exis- 
tence of Dirac points or the gapless phase is that \aj \ can 
form a triangle, or \ai\ + \aj\ ^ |afc|,{i,j,fc} = {1,2,3}. 
The conditions for the band collapse/dynamical localiza- 
tion become — aj = 0, S {1,2,3}. These condi- 
tions are two equations with two variables, which can be 
solved in general situation. 

Now we give a simple picture of the phenomena we 
have found. For an electric filed E(t) in arbitrary di- 
rection, the projection of this field in three directions 
along the bond (i.e., directions Ui,i — 1,2,3 (see figure 
1)) is Ei. Then the effective hopping constants along 
the bonds are V^^^ — V Ja{eEia/uj). Thus we may 
map our system to a system with honeycomb lattice and 
effective hopping constants V^^^ along bond u^. The dy- 
namical localization phenomena we have found are quite 
clear in this picture. In fact, for the case E//F and 
eEoa/uj = 4.81, = V^^^ = 0, the particle initially 

staying at site c can only oscillate between site c and site f 
(see Fig. 1(b)). It is the dynamical localization behavior 
shown in Fig. 4(a). For the situation eEoa/oj — 2.405, 
V^'^'^ — 0, the particle can only move through the path 
"abcde" shown in Fig. 1(b), which is the dynamical lo- 
calization in the Y direction show in Fig. 4(b). The 
occurrence of band gap and effective mass is also a conse- 
quence of the "bond ordering" , the dimerization pattern 
(see Fig. 1(b)). In the continuum limit, the fluctuation 
of hopping constant maps to axial vector potential [l5| . 
It is different from the honeycomb lattice with "Kekule 
pattern" , in which the fluctuation of hopping constant 
maps to the complex- valued Higgs field in the continuum 
limit. In that case the opening of energy gap is related to 
the spontaneous breaking of an effective axial U(l) sym- 
metry [l3| , and was proposed for the split of zero-energy 
Landau levels in graphene in the presence of magnetic 
field [l3l- By tuning the external electric field, one may 
effectively obtain different bond patterns. 

Realization in optical lattices The honeycomb op- 
tical lattice can be realized in experiment [J] by using 
three pairs of laser beams with co-planar propagating 
wavevectors ±qi, z = 1,2,3. The magnitude of the three 
pairs of wavevectors are the same and their directions 
form the angle of 120 degree with each other. The opti- 
cal potential has the form 

U = Y,U^ = UoY.cosip,■r), (6) 

i i 

where pi = qa - qa, P2 = qs - qi, and pg = qi - qa. 
Uq depends on the detuning between laser frequency and 
atomic transition frequency, and the resonant Rabi fre- 
quency proportional to the the quare root of laser in- 
tensity [l8[. Thus Uq can be tuned quite easily in ex- 
periments. For sufficient large detuning, we may neglect 
spontaneous scattering and the high frequency condition 



C/q/w << 1 can also be easily met. To simulate the 
electronic force in electronic system, we may let the co- 
ordinate of the mirror in each direction pi , reflecting the 
incoming traveling light wave, oscillate around its average 
with di cos{ujt). Its effect is to replace the atomic effec- 
tive potential Uq cos(pi -r) with Uq cos(pi -r — cos(wt)). 
In the moving frame of the potential, the atoms feel an 
effective force proportional to di cos{ujt), which is propor- 
tional to the electrical field component in the p^ direction. 

In summary, we have studied the quasienergy spectra 
and the related dynamics of a particle in honeycomb lat- 
tices driven by an AC field. By tuning the amplitude, fre- 
quency, and direction of the electric field, we can obtain 
quite rich phases, including gapped and gapless phases. 
Moreover we may tune the velocity of the "light" , and 
realize band collapse and the dimension reduction from 
2D to ID and OD (dynamical localization in one direc- 
tion and two directions). In the current paper, we have 
focused on the high frequency regime, i.e. V/lo << 1. In 
graphene with V around 2.8eV and the lattice constant 
a 1.4 A, our predictions may not be verified easily. How- 
ever, they should be quite easy to be observed in optical 
lattices. We leave the low frequency regime for future 
study. 
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